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We report that Bi2Se3 thin films can be epitaxially grown on SrTiC>3 substrates, which allow for 
very large tunablity in carrier density with a back-gate. The observed low field magnetoconductivity 
due to weak anti-localization (WAL) has a very weak gate-voltage dependence unless the electron 
density is reduced to very low values. Such a transition in WAL is correlated with unusual changes 
in longitudinal and Hall resistivities. Our results suggest much suppressed bulk conductivity at 
large negative gate- voltages and a possible role of surface states in the WAL phenomena. This work 
may pave a way for realizing three-dimensional topological insulators at ambient conditions. 
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A new state of matter, coined as topological insulator 
1 , , was recently discovered to exist in three dimensions 
3|, |4| . It is characterized by an insulating bulk and con- 
ducting surface states of massless helical Dirac fermions 
as a consequence of strong spin-orbit coupling (B"S| ■ The 
surface of the topological insulator can be viewed as a 
novel type of two-dimensional electron system (2DES) 
that may be a fertile ground for exploring exciting new 
physics as well as possible applications in spintronics and 
quantum information [2j|. Many of them are not possi- 
ble in its topologically trivial predecessors such as 2DES 
based on semiconductor heterostructures and graphene. 
Among all known materials that could host the topo- 
logical surface states, bismuth selenide (E^Ses) is most 
attractive because of a large band gap (~ 0.3 eV) in the 
bulk and a single Dirac cone in the surface energy spec- 
trum 

Despite numerous theoretical proposals that attempt 
to exploit their exotic properties, the 3D topological in- 
sulators have so far only confirmed to exist in ultra-high 
vacuum (UHV) environment 0, flfil - 14 1 . All of previously 
reported transport studies have been bothered by a con- 
ducting bulk [15rll8| . Even under circumstances of UHV, 
suppression of the bulk conductivity appears to be a chal- 
lenging task. Hsieh et al. reported that the chemical po- 
tential lies in the bulk band gap for a freshly cleaved, 
Ca-doped Bi 2 Se3 crystal, and an 18-hour stay in UHV 
would lift the chemical potential into the bulk conduc- 
tion band [l9| . Analytis et al. observed that an exposure 
of a cleaved Bi2Se3 surface in the air for 10 seconds is suf- 
ficient to introduce a significant amount of electrons into 
the bulk [l8(. It was demonstrated that the surface dop- 
ing effect can be overcome by tuning chemical potential 
with NO2 molecules deposited onto the Bi2Se3 surface in 



UHV so that a true topological insulator was obtained 
[l9| . This method is however cumbersome, if not impos- 
sible, for preparing devices suitable for transport studies. 
Therefore, it is very important to find more convenient 
means to lower the chemical potential in order to bring 
the surface states into the so-called topological transport 
regime, in which suppression of the bulk conductivity is 
a prerequisite. 

Here we report that large tunability in the chemical po- 
tential in Bi2Se3 is achieved with a back-gate. The gate- 
voltage dependencies of longitudinal and Hall resistivities 
as well as the weak anti-localization part of magnetocon- 
ductivity coherently point to a much depleted bulk and 
possible relevance of the surface states. 

The Bi2Se3 films were grown on SrTi03(lll) sub- 
strates with molecular beam epitaxy. Fig. la shows an 
in situ scanning tunneling microscopy (STM) image of 
a Bi2Se3 thin film. It displays many triangular-shaped 
flat terraces with a step height close to 1 nm, correspond- 
ing to one quintuple-layer (QL). The QL-by-QL growth 
mode of Bi 2 Se3 [20| enabled us to precisely control the 
film thickness. A 50 nm thick Pd film was subsequently 
deposited onto the back of the substrate to serve as a 
back-gate. The Bi2Se3 films were then patterned into 
Hall bars of millimeter sizes. The electron transport mea- 
surements were carried out in a cryostat with B up to 7 T 
and at temperature T — 1.8 K, unless otherwise specified. 
Hall measurements showed that undoped Bi2Se3 films are 
of n-type. We also found that doping the films with Ca 
can lower the electron density, consistent with previous 
work on doped bulk materials 
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The data to be pre- 
sented here were mainly taken from two samples. Sample 
A is a 48 nm thick undoped Bi2Se3 film, whereas sample 
B is doped with Ca and 10 nm thick. 

As shown in Fig. lb, the longitudinal resistivities at 
B = (p°x) °f samples A and B have similar tempera- 
ture dependence. When a perpendicular magnetic field 
is applied, both samples show positive magnetoresistance 
(MR) with a cusp-like minimum at B = 0. The MR of 
sample A in large fields has a parabolic independence, 
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FIG. 1: (color online) (a) STM image and low energy electron 
diffraction pattern (inset) of a Bi2Se3 single crystalline thin 
film grown on a SrTiC>3(lll) substrate; (b) T-dependence 
of p X x of samples A (lower) and B (upper) at B = 0; (c) 
Normalized magnetoresistance of samples A (lower) and B 
(upper); (d) Relative changes in carrier density with respect 
to Vg = extracted from Hall measurements (circles) and 
charging effects (squares). 



which is not observed in sample B. The Hall resistivities 
(pxy) of both samples varies linearly with B without any 
noticeable curvature. The electron densities extracted 



from the Hall coefficients are 

i!3 



i? a11 = -l/(eR H ) « 3.7 
and 1.1 xl0 13 cm~ 2 for samples A and B, respectively. 

The electron density can be modified substantially by 
applying the gate voltage Vg- Fig. Id shows the relative 
changes in the carrier density in sample B with respect 
to n s at Vq — 0. They were measured with a charging 
experiment, in which a constant current of 1 nA was used 
to charge the Bi 2 Se3/STO/Pd capacitor. The recorded 
charging time during the change of Vg thus provided a 
direct measurement of the variation in the total carrier 
density in the Bi2Se3 electrode. It has a nonlinear de- 
pendence on Vg, which is a typical dielectric response 
of STO single crystals to the applied electric fields 22 1. 
The relative change in the total carrier density with re- 
spect to n s at Vg = 0, An s , reaches ±2.5 x 10 13 cm~ 2 
at Vg — ±150 V, demonstrating a superb capability in 
tuning the carrier density. 

As shown in Fig. Id, the relative changes in n s ex- 
tracted from the Hall measurements, An^ a11 , are signif- 
icantly smaller in magnitude than An s . This suggests 
that the Hall measurements give underestimated n s val- 
ues, even though p xy has a perfectly linear dependence on 
B. This is not surprising because multiple types of charge 
carriers are very likely to coexist in the BioSgQ film. Ac- 
cording to the photoemission studies [Io|, 14, [3, fljj ]. 
the energy difference between the conduction band min- 
imum and the Dirac point of the surface states is only 
about 0.2 eV. This means that the two surfaces in the 
topological regime can only hold up to about 0.5 x 10 13 
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FIG. 2: (color online) Plotted as a function of An s are (a) 
p xx (hexagons) and its derivative (diamonds); (b) nf a11 (pen- 
tagons) and its derivative (triangles) . The negative gate volt- 
ages are partitioned into three zones (I, II, III) according to 
the degree of carrier depletion, (c) Calculated band diagrams 
based on a simple model described in the text. 



electrons/cm 2 . Therefore, even in sample B, which has a 
very low density, there exist a large number of conduction 
band electrons at Vg = 0. Considering that the bottom 
and top surfaces have different interfaces, we expect at 
least three different types of charge carries in the device, 
namely the bulk electrons and the carriers related to the 
top and the bottom surface states. Under this circum- 
stance, the Hall coefficients at the low B limit (lu c t <C 1, 
which is satisfied in both samples) do not provide accu- 
rate measurement of the total carrier density. It rather 
gives values in which more mobile charge carries have 
more contribution. 

Fig. 2 shows the carrier density dependence of longi- 
tudinal resistivity p xx an d rfjf . When n s is lowered, 
p xx increases monotonically, but its derivative, dp xx /e?n s , 
has a minimum at An s » —1.9 x 10 13 cm -2 , where the 
slope of dp^ x /dn s changes its sign from positive to nega- 
tive abruptly. Unusual features are also present in n^ a11 . 
The derivative dn^ a11 /dn s does not vary smoothly with 
n s . It rather shows sharp transitions at about —1.0 and 
— 1.9 x 10 13 cm -2 . The latter coincides with the transi- 
tion in dp^ x /dn s . 

Based on these observations, we divide the negative 
An s region in Fig. 2 into three zones to represent dif- 
ferent degrees of carrier depletion. In zone I, which 
spans from to —1.0 x 10 13 cm" 2 , p xx varies very slowly. 
The magnitude of dp^ x /dn s increases rapidly in zone II 
(—1.9 < An s < —1.0 xl0 13 cm~ 2 ). This is normally 
what one would expect from depleting a typical semi- 
conductor layer. At sufficiently low densities, metal- 
insulator-transition should take place in the bulk and the 
rapid increase in \dp1 x /dn s \ in zone II may be a precur- 
sor for such a transition. When An s < —1.9 x 10 13 cm -2 
(zone III), further depletion leads to the suppression of 
the bulk conductivity and hence the transport is dom- 
inated by the surface (or interface) states. This could 
account for the sudden change in the slope of dp^. K /dn s . 



Another noteworthy feature in this deep depletion regime 
(zone III) is that nf al1 only varies from -0.31 to -0.24 
xl0 13 cm~ 2 . Its slope becomes nearly zero when An s 
approaches to the lower limit. This would be difficult to 
explain if every part of the sample still remained n-type. 
A likely scenario is that the chemical potential at the 
bottom surface has dropped below the Dirac point. 

The profile of the chemical potential for a given An s 
can be calculated by solving the ID Poisson equation. 
For the simplest case, we do not consider band bend- 
ing near the top and bottom surfaces at Vg = 0. If 
we assume a bulk electron density of 1 x 10 19 cm~ 3 , 
Au s ~ —2 x 10 13 cm~ 2 is then required to reach the 
situation in which the chemical potential on the top sur- 
face is lowered to the conduction band minimum. Fig. 2c 
displays three band diagrams calculated for An s = —0.5, 
— 1.6 and —2.1 x 10 13 cm~ 2 , which may provide a crude 
picture for the evolution of the chemical potential and 
the carrier depletion as Vq decreases. Many experimen- 
tal factors, for instance, the band bending effects near 
the top and bottom surfaces as well as impurity levels, 
can certainly modify the results. Nevertheless, the basic 
trend shown in Fig. 2c should remain intact. 

Fig. 3 shows a detailed gate voltage dependence of mag- 
neto conductivity of samples A and B. It is defined as 
Aa(B) = cr xx (-B) — o"xx|s=o- F°r both samples Aa(B) 
has a cusp-like maximum at B — 0, which is typical for 
the weak anti- localization. For sample A (Fig. 3a), within 
the range of all accessible carrier densities, the low field 
part of Aa(B) varies very little. At B > 1 T, Aa(B) 
decreases as Vg becomes smaller. The high field part 
has a parabolic shape, which probably originates from 
semiclassical effects. 

The gate voltage dependence of the magnetoconductiv- 
ity of sample B is shown in Figs. 3b and 3c. For positive 
gate voltages, all Aa(B)-B curves are nearly identical, 
even though there are more than a factor of 2 change in 
Rh and an increase of 2.5 x 10 13 cm~ 2 in n s . In contrast, 
when Vg is lowered to —60 V and below, Aa(B) varies 
drastically with Vg ■ Measurements on several other sam- 
ples with various carrier densities that spread between 
those of samples A and B further confirm the following 
observation: the back-gate voltage can induce significant 
amount of change in the low field magnetoconductivity 
only when the carrier density is extremely low. 

In case of coexistence of a bulk Fermi surface with the 
surface states, the Bi2Se3 thin films may be treated as 
traditional 2D electron systems with strong spin-orbit in- 
teraction. Hikami-Larkin-Nagaoka (HLN) equation [23j 
is often applied to account for the magnetoconductiv- 
ity in this kind of systems. In the limit of very strong 
spin-orbit interaction and low mobility, i.e. <C t so and 
T0 <C t c , the HLN equation is reduced to 
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FIG. 3: (color online) Magnetoconductivity, defined as 
Aa(B) — a xx (B) — ct xx | s=0 of (a) sample A; (b) sample 
B at Vg =0, 20, 40, 60, 100, and 150 V; (c) sample B at 
V G = -150 V, -140, -130, -120, -110, -100, -80, and 
—60 V (from top to bottom). Fitted curves (lines) with Eq. (1) 
are also shown in (d) Vg = -150 V, (e) -110 V, and (f) 150 V 
for sample B. 



where r so (r e ) is the spin-orbit (elastic) scattering time, 
a = 1/2, ip is the digamma function, and B^ — 
h/ (ADer^) is a characteristic field related to the dephas- 
ing time r^. Here D is the diffusion constant. For the 
surface states in the topological regime, to the best of 
our knowledge, no theory has yet been given to quantita- 
tively describe the weak anti-localization arising from the 
Berry's phase of the Dirac fermions [1, 24|. We notice, 
however, that McCann et al. [25j developed a theory for 
graphene, in which WAL could be suppressed by inter- 
valley scatterings as well as chirality-breaking intravalley 
scatterings by lattice defects (e.g. ripples and disloca- 
tions). These effects are absent in the single-cone spin- 
helical system encountered here, so we expect that the 
magnetoconductivity of one topological surface of Bi2Sc3 
follows the same form as Eq. (1) with a — 1/2 [28]. Our 
fitting results with various samples showed that the con- 
ditions -C r so , t c are indeed satisfied at T — 1.8 K. As 
shown in Figs. 3d-f, Eq. (1) does provide reasonably good 
fits to the data at all gate voltages. 

Fig. 4a displays detailed fitting results for sample B. B^ 
varies very little for a wide range of densities, i.e. n s = 
— 1.4 to +2.5 x 10 13 cm 2 . When n s becomes lower (in zone 
II), Brf, starts to vary appreciably with Vq- It is striking 
that B^ nearly triples in a very narrow range of carrier 
densities in zone III. The correlation between the WAL 
data and p xx is shown more clearly in the B ( f,a xyi -An s 
plot (Fig. 4b), which has a minimum at the boundary 
between zones II and III. 
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FIG. 4: (color online) (a) (solid squares) and a (open 
circles), (b) B^a° x , as a function of An s - Here cr° x is the con- 
ductivity at B — 0. The partition of the negative Vg region 
is same as that in Fig. 2. (c) Calculated chemical potential on 
the top surface with respect to the conduction band minimum 
(CBM). The Dirac point (DP) is located at -0.2 eV. 



As shown in Fig. 4a, the coefficient a extracted with 
Eq. (1) is close to 1/2 for all carrier densities. One possi- 
ble explanation is that the Bi2Se3 thin film can be treated 
as a 2D system described by the HLN equation, but it is 
difficult to reconcile with the fact that increases by 
less than 15% for a wide range of densities (An s = — 1.4 



to 2.5 x 10 13 cm 2 ), while p xx and 



,HaIl 



are changed by 



a factor of 3.4 and 4.7, respectively. Such large changes 
are usually reflected in D and r^, and would change 
accordingly. 

An alternative explanation might come from treating 
the top and bottom surface states as well as the bulk elec- 
trons separately. Since both surfaces were measured si- 
multaneously in the experiment, this would lead to a = 1 
if they contribute equally. Apparently this is not the 
case. Considering that Vg is applied from the back of 
the sample, the electron states at and near the bottom 
surface will be influenced before any impact can be ex- 
erted on the top surface (see Fig. 2c). The weak carrier 
density dependence of B^ for such a wide range of Vg 
implies that the bottom surface or the nearby bulk does 
not contribute much to the magnetoconductivity at low 
fields. It rather mainly originates from the electron states 
at or near the top surface. This may be understood as 
follows. The bottom surface has many defects because 
of the lattice mismatch with the STO substrate. Due to 
the layered nature of the Bi 2 Se3 crystal and lattice relax- 
ation during the growth, the top surface is expected to 
be of considerably better quality. Consequently, the top 
surface has much higher carrier mobility (fi) and thus 
smaller B^ since D oc /z. This could result in negligi- 
ble contribution from the bottom surface. In addition, 
the Bi2Se3 crystal has inversion symmetry, which may 
account for insignificant contribution from the bulk. 

If the low field magnetoconductivity is indeed mainly 
contributed by the states at or near the top surface, the 
large increase B^ in zone III can be explained qualita- 
tively. In this low density regime, we anticipate a sig- 



nificant drop in the local chemical potential at the top 
surface as Vg becomes more negative (Fig. 4c). This will 
result in substantial decrease in the carrier density of the 
top surface states. For a 2D system the diffusion con- 
stant is D = l/2v pT e , where the Fermi velocity vp is 
constant for Dirac fermions. The larger B^ at lower den- 
sities can therefore be attributed to reduced screening 
and stronger electron-electron interaction effects, which 
are known to shorten r e and r^. In case that the con- 
duction band electrons are not fully depleted near the 
top surface, similar arguments can also be applied. The 
corresponding ^-dependence of B^ should, however, be 
qualitatively different because the dimensionality of the 
electron system is changed from 2D to 3D (or quasi-2D). 
The B^a^y. minimum observed near the boundary be- 
tween zones II-III in Fig. 4b might be associated with a 
transition of the electron states from a bulk-like to a more 
surface-dominating type. 

Finally, we briefly remark on the coefficient a, which 
has a sizable, density-dependent deviation from 1/2. 
This is not expected from the theories mentioned above. 
Possible explanations may arise from electron-electron in- 
teraction effects [26[ as well as states related to defects in 
Bi2Se3 jlo|. Furthermore, the large electric field induced 
by the back-gate voltage might significantly modify the 
spin-orbit interaction and hence spin relaxation via the 
Rashba effects [LJ, |27j . A thorough explanation of the 
data is out of the scope of this paper, but it certainly 
deserves further work along this direction. It will be of 
great interest to see whether such studies could provide 
a new method to distinguish a topological metal from an 
ordinary 2D electron system. 
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